o 

(N 
(N 



Mem. SA.It. Vol. 

©SAIt 2010 Memorie deiia * 




e ^OegliS9* 



Multiple stellar populations in the Globular 
Clusters NGC1851 and NGC6656 (M22). 

A. P. Milone^G. Piotto^L. R. Bedin 2 ,A. F. Marino 3 , Y. Momany 4 , and S. 

Villanova 5 

1 Dipartimento di Astronomia, Universita di Padova, Vicolo dell'Osservatorio 3, Padova, 
1-35122, Italy 

Space Telescope Science Institute, 3700 San Martin Drive, Baltimore, MD 21218, USA 
Max Plank Institute for Astrophysics, Postfach 1317, 85741, Garching, Germany 
European Southern Observatory, Alonso de Cordova 3107, Vitacura, Santiago, Chile 
Grupo de Astronomia, Departamento de Fsica, Casilla 160, Universidad de Conception, 
Chile 



Abstract. In the last years, photometric and spectroscopic evidence has demonstrated 

i-Q ■ that many, maybe all the Globular Clusters (GC) host multiple stellar populations. High- 

ft' resolution spectroscopy has established that, while most GCs are mono-metallic with no 

C*. [ significant abundance spread in j-elements, in all the globulars studied to date the presence 

5^ , of different stellar generation is inferred by the Na-O and the C-N anticorrelations. 

In this context, NGC 1851 and NGC 6656 are among the most intriguing clusters. Contrary 
to the majority of GCs, they host two groups of stars with different s-elements abundance 
that are clearly associated to the two distinct sub-giant and red-giant branches detected 
in their color-magnitude diagrams (CMD). In the case of NGC 6656 s-rich stars are also 
enriched in iron and calcium. Each s-element group exhibits its own Na-O and C-N anti- 
correlations thus indicating the presence of sub-populations and suggesting that the parent 
>0 , clusters have experienced a very complex star-formation history. In this paper we summa- 

rize the properties of multiple populations in NGC 1851 and NGC 6656. 

Key words. Stars: abundances - Stars: atmospheres - Stars: Population II - Galaxy: glob- 
r — . ular clusters 
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1. Introduction Spectroscopic studies have demonstrated 

that most GCs are mono-metallic with no de- 
In recent years, an increasing number of pho- tectable spread in their iron content and also 
/\ • tometric and spectroscopic observational evi- s-process elements do not exhibit large star- 

dence have shattered the paradigm of globulars to-star variations in the majority of globulars 
as the prototype of single, simple stellar popu- (e.g. Carretta et al. 2009a, D'Orazi et al. 2010 
lations (see Piotto 2009 for a recent review). and references therein). On the contrary, in all 

the clusters studied to date with large stellar 
samples, have been detected star-to-star vari- 
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Al (e.g. Carretta et al. 2009b, Pancino et al. 
2010 and references therein). These variations 
are related to correlations and anticorrelations, 
which indicate the occurrence of high tem- 
perature hydrogen-burning processes which in- 
clude CNO, NeNa, MgAl cycles and cannot 
occur in presently observed low mass GC stars. 

Today it is widely accepted that the ob- 
served light-elements variations provide strong 
support to the presence of multiple stellar pop- 
ulation in GCs with the second generations 
formed from the material polluted by a first 
generation of stars. On the contrary the debate 
on the nature of possible polluters is still open 
(e.g. D'Antona & Caloi al. 2004, Decressin et 
al. 2007). 

While abundance variations are well 
known since the early sixties, it was only the 
recent spectacular discovery of multiple se- 
quences in the CMD of several GCs that pro- 
vide an un-controversial prove of the pres- 
ence of multiple stellar populations in GCs and 
brought new interest and excitement in GCs 
research (e.g. Piotto et al. 2007). Photometric 
clues, often easy to detect simply by inspection 
of high-accuracy CMDs, arise in form of mul- 
tiple main sequences (MS, Bedin et al. 2004, 
Piotto et al. 2007, Milone et al. 2010), split 
sub-giant branch (SGB, Milone et al. 2008, 
Anderson et al. 2009, Piotto 2009), and multi- 
ple red-giant branch (RGB, Marino et al. 2008, 
Yong et al. 2008, Lee et al. 2009). 

Apparently the main properties of stellar 
generations, like the relative fraction of stars, 
their location in the CMD, the spatial distri- 
bution, and the chemical behavior differ from 
cluster to cluster. However, it is clear that some 
groups of clusters share all similar properties: 

1 . Na-O anticorrelation and C-N correlations 
are presents in all the GCs studied to date. 
When stars with available Na and O abun- 
dances have been identified in the U vs. 
(U-B) CMD, it was found that the group 
of Na-poor (O-rich) stars are spread on the 
blue side of the RGB, while the Na-rich (O- 
poor) population define a narrow sequence 
on the red RGB (e.g. Marino et al. 2008). 

2. In few cases the MS morphology supports 
the presence of stellar populations with dif- 



ferent helium. NGC 2808 has three dis- 
tinct MSs (Piotto et al. 2007) possibly as- 
sociated to three stellar population with 
primordial helium and with Y~0.33 and 
Y-0.40 (D'Antona et al. 2005). A spread 
MS have been detected in NGC 104 (47 
Tuc, Anderson et al. 2009) and NGC 6752 
(Milone et al. 2010) where there is also 
some hint of a MS split. 

3. In some GCs, like NGC 104, NGC 1851, 
NGC 5286, NGC 6388, and NGC 6656 
there is a double SGB (Milone et al. 2008, 
Piotto 2008, Anderson et al. 2009) associ- 
ated to two stellar groups with either a dif- 
ference in age by ~1 Gyr or with almost 
the same age but a significant difference in 
their overall C+N+O content (Cassisi et al. 
2008, Ventura et al. 2009, Di Criscienzo et 
al. 2010). 

4. Finally there is the 'extreme' case of to 
Centauri with either multiple MSs (e.g. 
Anderson 1997, Bedin et al. 2004), mul- 
tiple SGBs (e.g. Sollima et al. 2005, 
Villanova et al. 2007), multiple RGBs (Lee 
et al. 1999, Pancino et al. 2000) and large 
star-to-star variation in iron and s-elements 
(e.g. Johnson & Pilachowski 2010, Marino 
et al. 2010). Interestingly multiple popu- 
lations have been detected also in the Sgr 
dwarf galaxy central cluster NGC 6715 (e. 
g. Siegel et al. 2007, Piotto 2009 and ref- 
erence therein) and in Terzan 5 that is con- 
sidered the surviving remnant of one of the 
primordial building blocks that are thought 
to merge and form galaxy bulges (Ferraro 
etal. 2009). 

NGC 1851 and NGC 6656 are among the 
most studied clusters of the third group and are 
the main subject of this paper. In the following 
we will compare the observed CMDs of these 
GCs and the abundance of the critical elements 
Fe, O, Na, and s-elements. 



2. The color-magnitude diagram 

High-accuracy photometry obtained with the 
WFPC2 and WFC/ACS camera on board of the 
Hubble Space Telescope (HST) has revealed 
that the SGB of NGC 1851 is clearly split into 
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Fig. 1. Zoom of the CMD around the SGB for NGC 

6656 {upper panel) and NGC 1851 {lower panel) 
from WFC3/UVIS and ACS/WFC data (GO10775, 
P. I. Sarajedini and GQ12311, P. I. Piotto). 



two branches with the bright SGB component 
containing about the 60% of the total num- 
ber of SGB stars and the remaining 40% of 
stars belonging to the faint SGB (Milone et al. 
2008). Similarly to the case of NGC 1851, af- 
ter a careful correction for differential redden- 
ing, the SGB of NGC 6656 also revealed a bi- 
modality with a fraction of bright over faint- 
SGB stars similar to the one measured in NGC 
1851 (Piotto 2009, Marino et al. 2009). A 
comparison of the SGBs of NGC 1851 and 
NGC 6656 is given in Fig. Q] where we show 
a zoom of the WFC3/UVIS CMDs around the 
SGB. The spread SGB has been explained in 
terms of two stellar generations, only slightly 
differing in age, the younger one having an 
increased C+N+O abundance (Cassisi et al. 
2008, Ventura et al. 2009). As an alternative 
possibility, we note that observations are con- 
sistent with two stellar populations with con- 
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Fig. 2. CMD of NGC 6656 from ground-based pho- 
tometry corrected for differential reddening. 



stant CNO but differing in age by ~1 Gyr 
(Milone et al. 2008). Significant star-to-star 
variations in the overall CNO abundance have 
been detected among NGC 6656 RGB stars by 
Marino et al. (2010) and in two out four NGC 
1851 giants by Yong et al. (2009). The letter 
results is not confirmed by the recent analysis 
by Villanova et al. (2010) who did not find sig- 
nificant CNO abundance variation in 15 RGB 
stars of NGC 1851. 

Neither the MS nor the RGB of these clus- 
ters shows any hint of bimodality or spread in 
the wp606w -fn-pnAw color. Ritcher et al. (1999) 
found that NGC 6566 exhibits a bimodal distri- 
bution in the m\ index among RGB stars and a 
similar bimodality was detected by Calamida 
et al. (2007) for the case of NGC 1 85 1 . A sim- 
ilar split among RGB stars of both clusters has 
been observed in the hk index by Lee et al. 
(2009) and in the U vs. (U-I) and Uvs. (U-V) 
CMD by Han et al. (2009) and Momany et al. 
(2004, Fig. |2]). In these cases the RGB compo- 
nents are clearly associated to the two SGBs. 



3. Chemical abundances 

In this section we investigate the behavior 
of proton-capture elements, iron, sodium and 
oxygen abundances as they provide crucial 
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Fig. 3. [Y/Fe] (upper panels) and [Fe/H] (lower 
panels) as a function of [Ba/Fe]. Filled circles and 
empty triangles represent the .s-rich and the s-poor 
stars. The red crosses with error bars indicate the 
average abundance of stars in each group. 



information on multiple stellar populations. 
A peculiar property of the cluster pair NGC 
1851-NGC 6656 is the large scatter in the 
abundance of those n-capture elements that are 
associated to ^-processes. This result comes 
from Marino et al. (2009, 2010) for NGC 6656 
and Hesser et al. (1982), Yong et al. (2008), 
and Villanova et al. (2010) for NGC 1851. As 
an example, Fig. [3] shows the abundance of 
[Ba/Fe] as a function of [Y/Fe] for NGC 6656 
(Left Panel) and NGC 1851 (Right Panel). In 
both clusters n-capture elements are clearly 
segregated around two distinct values of bar- 
ium and yttrium in sharp contrast from what 
found in most GCs where the abundance of 
these elements does not exhibit significant star- 
to-star variations. 

Since the solar-system abundances of bar- 
ium and yttrium are mainly due to s-process, in 
the following we will call stars with high and 
low [Ba/Fe] ([Y/Fe]) as 's-rich' and 's-poor'. 
Following Marino et al. (2009), we will also 
use different symbols in Fig. [3] and in the sub- 
sequent figures to separate these two groups of 
stars. 

When we compare the iron abundance for 
NGC 6656 with the ^-elements abundance, we 



find a strong correlation, with s-nch stars hav- 
ing systematically higher [Fe/H] as shown in 
the lower-left panel of Fig. [3] The s-nch and s- 
poor groups have average [Fe/H] r ; c h = -1.68 + 
0.02 and [Fe/H] poor = -1.82 + 0.02, respec- 
tively (5[Fe/H] = 0.14 + 0.03, Marino et al. 
2009). This result demonstrated that, at odds 
with 'normal' monometallic GCs the different 
stellar populations of NGC 6656 have signifi- 
cant difference in their iron content. 

The presence of intrinsic spread in iron in 
NGC 1851 has been widely discussed in liter- 
ature. Yong et al. (2008) analyzed UVES spec- 
tra for eight RGB stars and detected a dis- 
persion in [Fe/H] of ~0.10 dex suggesting a 
possible intrinsic spread in iron. A small but 
detectable metallicity spread, compatible with 
the presence of two groups of stars with a 
metallicity difference of ~0.07 dex has been 
claimed also by Carretta et al. 2010 on the basis 
of their abundance analysis of more than 120 
RGB stars. These results are not confirmed by 
Villanova et al. (2010) who studied high S/N 
MIKE/Magellan spectra for 15 RGB and found 
a smaller [Fe/H] spread compatible with obser- 
vational errors. Villanova et al. (2010) detected 
no significant difference in iron between the s- 
rich and s-poor as shown in the lower-left panel 
of Fig. S 

In the light of these results we matched 
spectroscopic data on NGC 6656 from Marino 
et al. (2009, 2010) with the Stromgren photom- 
etry by Richter et al. (1999). Results are shown 
in Fig.|4]where we plotted the m\ versus (b-y) 
diagram for NGC 6656. At odds with what ob- 
served by using the V vs (V-I) diagram, where 
the stars belonging to the two groups appear to 
populate an unique RGB sequence, it is evident 
that s-rich and s-poor stars define two distinct 
RGBs. As the mi index is strongly dependent 
by the CN bands strength, we expect this bi- 
modality as due to the overabundance in both 
C and N measured by Marino et al. (2010) in 
.s-rich stars. 

3.1. Na-0 anticorrelation 

Sodium versus oxygen is plotted in Fig. [5] for 
the pair NGC 1851-NGC 6656. The most im- 
portant result is that in these clusters both the 
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Fig. 4. the mi, (b-y) diagram for NGC 6656 cor- 
rected for differential reddening, .v-rich and j-poor 
stars are represented with full circles and open tri- 
angles. 



s-rich and the s-poor group have its own Na-O 
anticorrelation. In the case of NGC 6656 also 
the C-N anticorrelation has been detected in 
both the s-groups by Marino et al. (2010). 

In both clusters s-rich stars have, on aver- 
age, higher sodium but while in the NGC 6656 
both s-poor and s-rich stars span almost the 
same [O/Fe] range, in NGC 1851 s-rich stars 
are extended toward lower oxygen suggesting 
that in this case the interstellar medium could 
have been polluted by more massive stars. 

4. Conclusions 

We have summarized some spectroscopic and 
photometric findings on multiple stellar popu- 
lations in the GCs NGC 1851 and NGC 6656. 
From spectroscopic and photometric ob- 
servations, it comes out that the SGB and 
RGB split is related to the presence of two 
groups of stars with different ^-element con- 
tents and a possible difference in CNO abun- 
dance. According to this interpretation, s-poor 
stars would constitute the first population. The 
second stellar generation should have been 
formed after the AGB winds of this first stel- 
lar generation had polluted the protocluster in- 
terstellar medium with ^-elements. In the case 
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Fig. 5. Na-0 anticorrelation for NGC 6656 (upper- 
left panels) and NGC 1851 (upper-right panel) RGB 
stars, s-rich and .s-poor stars are represented as full 
circles and open triangles. In the upper panels the 
fiducial lines of the Na-0 anticorrelations for the s- 
rich and the .v-poor stars of NGC 6656 are as contin- 
uous and dashed gray lines respectively. In the lower 
panels we plotted [Na/H] as a function of [O/H]. 



of NGC 6656, this second generation may have 
formed from material that was also enriched by 
core-collapse supernovae ejecta, as indicated 
by their higher iron, magnesium, and silicon 
content, and the lack of correlation of the iron 
content with a pure r-process element (Eu). 

This scenario is seriously challenged by the 
recent findings on light elements Na, O, C, and 
N. As already mentioned in Sect. 1, the Na- 
O and the C-N anticorrelations can be consid- 
ered the signature of multiple star-formation 
episodes. This implies that both clusters are 
composed by two groups of stars with distinct 
s-element content that define double SGB and 
RGB observed in the CMD but each group is 
made by multiple sub-populations of stars with 
different Na and O abundance. NGC 1851 and 
NGC 6656 do not host only two stellar popula- 
tions but have experienced a much more com- 
plex star-formation history with similarities to 
the 'extreme' cases of a> Centauri (Da Costa et 
al. 2009, Da Costa & Marino, 2010) and the 
Sgr dwarf galaxy central cluster NGC 6715. A 
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detailed comparison of the abundance patterns 
in a> Centauri and NGC 6715 with those of 
NGC 1851 and NGC 6656 is strongly needed 
to resolve this problem. 
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